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1. Introduction 

The problem of finding optimal radiation doses for organs or tissues in therapy of cancer belongs to 
the principal ones in the modern Radiobiology (see, e.g., [20], US], HB|, GE]i [50] - anc ^ references 
therein). The complexity of this problem is related to nondeterministic character of the oncological 
therapy results. 

The aim of this paper is to study stochastic models for collective effects in the behaviour of the irradi- 
ated cells. We provide a generalization of the well-known critical volume (CV) normal tissue complication 
probability (NTCP) model to comprise a concept of functional subunits (FSUs) dependence. Such a model 
is beyond the scope of 9 and for its investigation new limit theorems are required. Note that here the 
links between Probability and Geometry are stipulated by the dependence structure of a random field 
under consideration which is governed by the configuration of a graph used as a parameter set. Moreover, 
it seems natural from the biological view-point to assume some dependence in the collective performance 
of cells (and FSUs). 

The description of the FSUs behaviour by means of non-binary random variables is considered as 
well. Other important biological response models and further research directions are tackled in the last 
Section. 



2. Accuracy of the CV NTCP model 

We recall the basic critical volume model (see ^S]> EH) an( i after that consider more carefully its 
framework. The organ (or tissue) modelled is assumed to be composed of independent FSUs and it 
is supposed that complications in its functioning arise only if sufficiently many FSUs ("the functional 
reserve") are destroyed. More precisely (see, e.g., j2Z])) the assessment of the impact of irradiation 
is divided into two stages. The first one is the reaction of the cells forming an FSU which gives the 
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probability of death or survival for the FSU. The second one is a probability to have sufficient number of 
FSUs survived the irradiation in order to maintain the organ's functionality (it means the complications 
absence probability). Note that the serial and tumor control models are special cases when the functional 
reserve consists of one or all the FSUs respectively. One speaks also of the critical element models. The 
last special models are natural for organs having the one-dimensional shape (e.g., spinal chord). For 
other type of organs it is desirable to consider the " parallel architecture" . 

For the sake of simplicity we start with a model of organ consisting of n similar FSUs. Introduce in- 
dependent identically distributed (i.i.d.) random variables (r.v's) X\, . . . , X n defined on some probability 
space (fi, !F, P) such that 

P(Xi = l)=p, P(Xi = 0) = q where q = 1 - p, < p < 1. (1) 

More exactly consider an array of r.v's X,- t = X,(D) where D is a positive parameter (the irradiation 
dose). Assume that Xi(D) shows the state of the z-th FSU after irradiation of dose D (i = 1, . . . ,n). 
Namely, the random event {X t = 1} means that the z-th FSU is killed and {Xi — 0} corresponds to the 
case that the i-th FSU survives. As usual the argument uo € ft is omitted and we write Xi(D) (or simply 
Xi) instead of Xi(u>,D). 

Set S n (D) — Y^i=i Xi(D). In other words consider a random variable equal to the number of FSUs 
killed due to irradiation dose D. Thus NTCP is P(S n (D) > L) with the threshold L being some positive 
integer. 

Using the convergence rate estimate in the central limit theorem (see, e.g., ^UJ; P- 323) one has 

sup \P{S n {D) > x) - 1 + $(z)| < c/y/np(D)q(D) (2) 

— oo<x<oo 

where z — (x — np(D)) / <J np(D)q(D) , p(D) and q{D) appear in Q for Xi = Xi(D), positive constant 
c < 0.7975 and the c.d.f. of a standard normal law 

1 f z _„i 
$(z) = -= / e 2 du, z S R. 

V27T J-oo 

Thus if np(D)q(D) is large enough, 2 for any threshold x (possibly depending on n and p(D)) the 
following relation holds 

P(S n (D) > x) w 1 - $(z) (3) 

where z is introduced above and the proximity is evaluated by the right hand side of . It is well known 
that the power — | of n in J5J is the best possible (see, e.g., [221 f° r that and also about non uniform 
estimates of the convergence rates in CLT for independent summands) . 

Consequently, for a given value 7 S (0,1) (close enough to 1), using the table of function $ one 
can find z = z 7 as the unique root of the equation $(z) = 7. Then the probability of complications is 
approximately 1 — 7 (with exactly specified boundaries) and the threshold 

x 1 = np(D) + {np{D)q{D)) l ' 2 z 1 . (4) 

Evidently, £ 7 = x 7 (n,p(D)). 

The distribution of r.v. S n has atoms at the points 0, . . . , n. Therefore one could choose an integer 
threshold L 1 = [x 7 ] where [•] stands for the integer part of a real number. In this case a trivial estimate 
for the continuity module of the function $ shows that 

|P(S„(D) > L,) - 1 + *W| < -^===p, (c + -±=) . 

We observe that if np(D)q(D) is large enough then a search for an integer threshold L 7 is not important. 

2 Otherwise one has to use different approximations. 
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Note in passing that for the number of " successes" S n in the Bernoulli scheme (of n independent trials 
with probability p for "success") one can apply (see, e.g., a little bit different approximation using 
for < k < m < n the relation 

P(fc < S n < m) = $(t 3 ) - + -^-{(1 - t 2 )e- t2/2 }l^ + A 

OV 27T(T 



where tr = ^npq, t\ = (fc — 5 — np)/a, ti = (m + 5 — np)/a and the error term A satisfies for cr > 5 the 
inequality 



|A| < (0.12 + 0.18|p- g |)cr- 2 + e 



-3<t/2 



Actually we deal with an equivalent description of the well-known critical volume model. Namely, 
suppose that the volume 3 of the irradiated organ is V and let Vt represent the volume of the i-th FSU 
(i = 1, . . . ,n). Clearly instead of S n (D) (the number of killed FSUs) we could consider the random 
damage volume V — Yli-x =1 ^ n ^ ne case when Vi — V/n, i = 1, . . . , n, one has 

V = V S -^1. (5) 

n 

Thus we come to description of the irradiation result in terms of the damage volume and one can specify 
the threshold v c for P(V > v c ). 

Formula suggests that it is natural to introduce a threshold of the type 4 x = xn where xr G (0, 1) 
is the fraction of killed FSUs. Thus 

x = p + c(p(l — p)) 1 ^ 2 where p = p(D), c = z 1 n~ 1 ^ 2 . (6) 

Note that z 7 > for 7 > 1/2 and consequently c > 0. Evidently for c = (i.e. 7 = 1/2) one has 
x = p. For each c > the graph of a function x = x(p) has the following features. One can easily verify 
that x(pi) = 1 for p\ = 1/(1 + c 2 ), x 1 (0+) = +00, x! {1— ) = —00 and the concave function x(p) attains 
its maximum te* = |(1 + \/l + c 2 ) at the point p* = |(1 + 5 ). Moreover, for 7 > 1/2 one has 

< inf (x(p) — p) < sup (x(p) — p) < z 1 j{2\fn). 

0<P<1 0<p<l 

Note also that if p(D) > p\ then relation © can not be satisfied for any x € (0, 1). 

On the other hand, given n G N, 7 > 1/2 (i.e. c = z^/^/n) and >f £ (0, 1), there is a unique root p = p 
of equation (JSJl 

P= fx+^ + c^->, 2 + ^ 7 V(! + c2 )- ( 7 ) 



We remark also that p(D) should be nondecreasing function on (0, 00) and if p(D) is continuous then 
for any x £ (xi, X2) where x% = x(ini o>o p(D)) and xi = min{l, x(sup D>0 p(D))} there exists (unique 
if p is strictly increasing) D such that p(D) = p. 

Now we discuss the models providing p{D) . Assume that every FSU consists of no cells. The surviving 
fraction of these cells after irradiation of dose D is determined (see, e.g., [HO]) by 

SF(D) = exp{-aD} 

where a > is the radiosensivity 5 of the cells. This is a so-called single-hit model. Suppose that each 
cell of FSU behaves in the same manner as other ones. Usually one admits that an FSU can regenerate 
from a single surviving cell, which means it is disabled only when no cell survives. Thus the probability 
of killing an FSU due to irradiation of dose D is 

p{D) = (l-e- QjD ) no . (8) 



3 Possibly a length or an area, the interpretation depends on the model of an organ 
4 We do not use in this paper the theory of large deviations for sums of r.v's. 
5 One writes also SF(D) = exp{— D/Dq} where Dq is called the mean lethal dose. 
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To obtain (jSJ) one supposes that all no cells in a FSU evolve independently of each other. 

Now assume that every cell contains m targets, each of them must be hit at least once to inactivate 
the cell. Then the probability that all targets of a cell will be hit at least once is (1 — e~ aD ) m . Thus for 
this multi-target model 

SF(D) = 1 - (1 - e - aD ) m , p{D) = (1 - e - aD ) vmo . 

Most experimental survival curves have an initial slope whereas the multi-target/single-hit model predicts 
no initial slope. To have a more adequate description one uses the family of functions 

S'F(D) = e- QD (l-(l-e- /3D ) m ), a>0, 0>O, 

Note (see, e.g., ^1) that for SF(D) a linear quadratic (LQ) model is also widely used with 

SF(D) = e - {aD+fm ' 2 \ a>0, /3>0. 

Now we concentrate on the generalization of the CV model considered above. 



3. Variant of the central limit theorem for dependent random fields 

Let X(D) = {Xj(D),j G Z d } (d > 1) be a family of random fields defined on a probability space 
(f2, T, P) for D > 0. Employing instead of an integer lattice Z d a parameter set T = <5Z d , with S > 0, 
one can easily reformulate all the results for a family of FSUs assuming, e.g., that the fc-th FSU is a cube 
with a center at a point tk G T and with an edge length equal to S. In other words one can use a scale 
appropriate to the problem under consideration. Thus without loss of generality we restrict ourselves to 
the study of a random field X(D) on a lattice 7L d . Moreover, we can assume that the random variable 
Xj(D) describes the state of the corresponding FSU after its irradiation of dose D. This gives us a 
possibility to consider not only the death and survival of an FSU but also to consider the " intermediate" 
states. Then the general (collective) effect of irradiation is represented by the following sum 

S(U,D) = Y,Xi(D) 
jeu 

where U is a finite subset of Z d . For a fixed D we also write simply Xj and S(U). 

For a finite set I C Z d with cardinality |/| introduce the a-algebra A(I) = a{Xj,j G 7}, that is 
consider the tr-field generated by a field X over a set I. 

There are different methods (see, e.g., 0, ^SJ) to describe the dependence structure of a field X. 
Here we use the maximum correlation coefficient for A(I) and A(J) over finite disjoint sets I, J C Z d 
which is defined as follows 

p(I, J) =sup{|c rr(£,77)| : { e L 2 (0,^(/), P), ?7 <= L 2 (£l, A(J), P)} (9) 

where corr(^, rf) is the correlation coefficient for (nondegenerate, square integrable) real-valued random 
variables £ and ij measurable with respect to cr-algebras A(I) and A(J). 

Assume that for all finite disjoint sets I, J C Z d , some positive Co and A one has 

p(I,J) <c \I\\J\(dt S t(I 1 J))- x (10) 

where 

dist(I, J) — min{||<2 — j\\ : q<El,j<E,J}, \\z\\ = max \Zk\, z € 

l<k<d 

Remark 1. Employing condition has the following motivation. In many stochastic models it 
is reasonable to assume that dependence between the random variables {Xj,j G 1} and {Xj,j G J} is 
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rather small if the distance between / and J is large enough. However, due to the paper |12j it was 
realized that for random fields (in contrast to stochastic processes corresponding to the case d = 1) 
one cannot, in general, measure the dependence between A(I) and A(J) only in terms of the distance 
between I and J. Namely in many situations the dependence between A(I) and A{ J) could increase for 
sets I and J growing, e.g., in such a way that the distance between them is fixed. Dependence notions 
based on correlations are quite familiar in various domains of applied probability. Appearance of the 
factors \U\, \V\ and parameters c, A in IjlUfl is intended to account, in a qualitative sense, for the effect 
of possible increase of dependence between er-algebras A(I) and A(J) when / and J are growing so that 
the dist(I, J) is preserved. Besides, (|10|) implies the same (i.e. power-type) decrease of correlations when 
I and J are moved apart so that the distance between them tends to infinity. Of course, a simple and 
natural hypothesis of m-dependence is a particular case of our condition. Recall that a random field 
X = {Xj,j 6 Z d } is m-dependent (with some to > 0) whenever A(I) and A(J) are independent if 
distil, J) > to. Thus, we include, in particular, a useful model of dependent nearest neighbouring FSUs. 
See also Remark 2 and Section 5. 

First of all we establish the central limit theorem (CLT) with convergence rate for partial sums 

S(Un) = X 3> ne N ' ( H ) 

of multi-indexed dependent r.v's where summation is carried over the integer cubes U n = [—n, n] d D Z d , 
n e N. 

Theorem 1 Let X{D) = {Xj(D),j S Z d }, D > 0, be a family of strictly stationary random fields such 
that for some S 6 (0, 1], C2+s(D) > and any D > 

E\X (D)f+ s <c 2+s (D). (12) 

Assume that condition (|1(J|) holds for all fields X(D) with the same A > 4d(l + S)/S and cq. Then there 
exists v = v(d, \,5) > such that for each D > and any neN 

sup|P((S(t^,D) - \U n \EX (D))/(o-(D)\U n \^ 2 ) < x) - $(x)| < A\U n \~ v (13) 

where A = A (d, A) max{l, coEX$(D)}max{l, c 2+S (D)/(T 2+s (D)} and 

a 2 {D) = cov(X (D),Xj(D)) ? 0. (14) 

Proof is based on the classical blocks technique initiated by Bernstein, so we only indicate the main 
steps and concentrate in the next Section on a statistical version of this result. 

For every n € N introduce p = p(n) = [n a ] and q = q(n) = [n 13 ] where < (3 < a < 1 and [•] 
stands for an integer part of a number. Consider k — kin) — [{In + l)/(2p + q)]. Then one can write 
\—n, n] = ii U I[ U . . . Ik U I' k U I' k ' where I m ,I' m , I' k ' are disjoint intervals of the form I m = [a m , a m + 2p], 
I' m = {a m + 2p, a m + 2p + q) and I' k ' = (a k + 2p + q,n] (I' k ' can be empty, a% = —n,a m G [— n, n], 
to = 1, . . . , k). 

Set Bi = I it x . . . x I id n 1 d where i — (ii, . . . , id) € M n = {1, . . . , k} d and let V n = Ui e M n Bi- 
It is easy to verify that for all n large enough 

^HS(U n ) - ES(U n ) ~ (S(V n ) - ES(V n ))) 2 < - g Ico^Xo,^)! < 4dn"^(D) (15) 

where 7 = min{l — a, a — (3} and the series 



\cov(X (D),X j (D))\=v(D) 



(16) 



converges in view of l|l(J|l for A > d. Consequently, 



S(U n )-\U n \EX S(V n ) - \V n \EX 



<r(D)\U n \V* 



< 2(dv{D)) 1 ^ 2 n-^ 2 a- 1 (D). 



(17) 



Lemma 1 Let X{D) = {Xj(D),j £ Z d } be a wide -sense stationary random field such that (|10f> holds 
with some A > d. Then for all n £ N 



varS(U n ) 



\U n 



-o*{D) 



where a = ao(d, X)cqEXq(D) and 



< af(n,d,X) 



d < X < d+ 1, 



(18) 



Proof. One has 



f(n,d,X) = { (l + hm)/n, X = d + 1, 

A > + 1. 

-5(17 W ) -<7 2 (£>)|tf n || < £ £ law^.^J^Ti + Ta 



where the sums T\ and T 2 are taken respectively over the sets {i £ U n ,j (f. U n , \\i — j\\} — r,r > n} and 
{i £ U n ,j U n , =r, r £ {1,.. .,n}}. Evidently, 

Ti < 2d3 d - 1 \U n \c Q EX 2 ^ r d - 1 - x < 2d3 d - 1 c EX 2 n d - x \U n \/ (X - d), 



T 2 <Y, E IcoviX^X^l 

r—1 n—r<\\i\\<n,\\j—i\\=r 
n n 

< 2dc EX 2 E r- x {{2n + l) d - (2(n - r) + l) d )(2r + l)*" 1 < 4d 2 3 d - 1 c (2n + l^EX 2 ^ r d_ 

r=l 

Using a trivial estimate 



r=l 



71 f-Tl 

Yr-^<1+ x'^dx, [i>0, 

r=l , ' 1 



we come to relation <|18f) . The Lemma is proved. 

Set Yj = S(Bj) — ES(Bj) where Bj belongs to a collection of "large" blocks, j £ M n = {1, . . . , k}, 
k = k(n). Clearly Yj = Yj(p n ,D), Introduce independent copies Zj, j £ M n , of random variables Yj, 
j £ M n . 

Then it is easily seen that for any i 6 R, i 2 = -1 and all n large enough 

|Eexp{zi J2 ^}-Eexp{z< £ ^}| < 4|Af„|c q- A (2p + l) d |ig < 4c 3VV\ 



Thus 



E exp 



MD)\Un\^ 2 

( itJ2jeM n Yj 
\a{D)\U n \V* 



<v(D)/a 2 (D), E 



< v{D)/a 2 (D). 



E exp 



<?{D)\Un\ 1 ' 2 , 

< mm{4c 3 d n 2d q- x ,2\t\y/ : 4l))/a(D)}. (19) 



G 



Using Lemma 1 one can verify that for all n large enough 



< {Adn-"> + aa- 2 (D)f(p, d, A)) 2 



\v{D)\U n \V* (wE iS M„^) V2 , 
where 7, a and / are the same as in (|15f) and (|18p. Therefore for A > d + 1 we can write 



^jeM n Zj 



*(D)\U n \V* (varE jeMn Zj) 1/2 

where Ci = Ci(d, a) max{l, <7~ 2 (_D)}, r = min{7, a}. 

Now the Esseen inequality implies that for every T > one has 



< c in - T 



(20) 



\t\<T 



Eexpj^ ^^^J - exp{-^} 



(it + a^T 



(21) 



where a\ and 02 are absolute positive constants. 

Applying the Berry-Esseen estimate of the convergence rate in the CLT for independent summands 
Zi,i £ M n , with finite absolute moments of order 2 + S (see, e.g., ^Ujj P- 322), using (|T7J) - (|2U)l and 
estimating the integral in the right hand side of i|21[) as a sum of integrals J| t | <:1 y T and J 1 / T< | t |< T and 
finally taking T — brfi with appropriately small £ and specified b > we arrive at This completes 

the proof of Theorem 1 . 

Remark 2. There are many versions of the CLT for random fields under various dependence condi- 
tions (see, e.g., 0, 0, 0, [TT] . 0) and references therein). In the same manner we could use instead 
of the maximal correlation coefficient p, e.g., the Rosenblatt-type mixing coefficient. We proved here the 
CLT with rate because it permits to establish the law of the iterated logarithm (announced in under 
the dependence conditions of the type (|10|l . It is worth mentioning that to this end we need only arbitrary 
slow power- type estimate of the convergence rate in the CLT without specifying an exponent v in 113|) . 
We do not provide here an explicit cumbersome expression for v. More restrictive mixing conditions than 
BlOll . i.e. p(I, J) < cq\I\\J\ exp{— a dist(I, J)} where a and b are some positive parameters, were recently 
used in [23 , Q5| ( see also the references therein) for CLT and LIL. We do not consider here growing 
subsets U n C Z d more general than " integer" cubes. For generalizations of this kind we refer to , [7] • 



4. Statistical version of the CLT 

There are two ways for applications of Theorem 1. Namely, if we believe in the model describing the 
stochastic behaviour of each FSU (see Section 1) then we can calculate EXq. However, the problem for 
dependent FSUs is the following one. Now we cannot claim (in general) that the variance of the sum 
S(U n , D) is equal to the sum of variances of summands. Thus for every D > 0, in contrast to the CLT 
for the Bernoulli scheme, i.e. 

S n (D)-np(D) v 

■ Z ~ N (0, 1J as n — > 00 



y/np(D)(l-p{D)) 

discussed in Section 1 (here p(D) = EXo, Z is a standard normal r.v.), the relation 

S(U n ,D)-\U n \EX (D) v 



a{D)\U n \^ 



iV(0,l) as rnoo (22) 



v 



contains an unknown function u{D). As usual — ► stands for weak convergence of random variables 
distributions. 
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Fortunately it is possible to construct a sequence of nonnegative statistical estimates C(U n ,D) for 
cr 2 (D) such that for any D > 

C(U n , D) ^ a 2 {D) as moo (23) 

p 

where — > means the convergence in probability as usual. We employ here a family of consistent statistical 
estimates introduced in |Sj for random fields 6 , for stochastic processes we refer to the paper |21j . 
Then by virtue of (|22|l and ll'-'il) we come, for every D > (if <j(D) 2 ^ 0), to the formula 

(C(U n ,D)\U n \)- 1 / 2 (S(U n ,D)-\U n \EX (D))°Z~N(0 1 l) as n -> oo. (24) 
In other words a random normalization is used in the CLT. 

Consequently to determine (approximately) for a given value 7 € (0, 1) the threshold x 7 we can apply 
the following analogue of formula J3j) 

P(S n >i)«l- *(0 (25) 



where £ = (x — nEX (D))/y C(U n , D)\U n \, x £ R, n G N. However, now in the right hand side of 
there is a r.v. $(£), i.e. we use $(£) as statisitical estimate for P(S n < x). Note that we have used only 
the value EA (-D) provided by the model of stochastic behaviour for FSUs and we did not suppose here 
that the collective effect of the evolution of cells under irradiation is described by independent binary 
random variables. 

Another way of using Theorem 1 is to construct approximate confidence intervals for the unknown 
mean value EXq(D) without hypotheses concerning the explicit formulas (discussed in Section 1) for 
distribution of random variables Xj, j G U n . 

Thus in both cases it is desirable to establish the CLT for dependent random fields using random 
normalization. 

Remark 3. As far as we know, in previous applications of the CLT to NTCP models for independent 
FSUs the question of convergence rate was not raised, so Section 1 covers this gap. However, the same 
question in case of dependent FSUs is more involved. We intend to investigate the accuracy of the 
proposed model in a special publication. One can consider Theorem 1 as the first step in this direction. 
Moreover, we can obtain the power-type estimate in the CLT with random normalization. However, the 
rate of convergence will be slower than that for independent random summands. The effect of convergence 
rate in the CLT sensitivity to the dependence conditions was demonstrated for positively or negatively 
associated random fields in 0]. 

For j G U C 7L d (1 < \V\ < 00) and b = b{U) > set 

Kj{b) ={teZ d :\\j- t\\ < 6}, Qj = Qj{U, b) = Un K^b), (26) 



Note that the averaged variables S(Qj,D)/\Qj\ arise for dependent summands (in contrast to the 
traditional estimates of variance used for independent observations) . 

Theorem 2 Let the conditions of Theorem^be satisfied. Let {J7 n } n eN be a sequence of "integer" cubes, 
i.e. U n — [— n,n] d D "L d , n G N (d > 1). Assume b(U n ) = b n where {fe„}„gN is a sequence of positive 
integers such that 

b n — > 00, b n — o(n) as n — > 00. (28) 
Then for every D > relation holds if a(D) ^ 0. 

6 Vector- valued random fields satisfying other dependence conditions are studied in 
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Proof. The estimate C(U,D) introduced by means of (|27|) and <J 2 {D) are invariant under the trans- 
formation Xj i !■ Xj — EXq, j G 7L d . So, without loss of generality we can further on assume that EXq = 0. 
Let stand for the norm of a real-valued random variable £ in a space L = L (£l,!F, P). For any 

U n C Z d and n€fi one has 



\\C(U n , D) - a 2 (D)\\ L < h(U n) D) + I 2 (U n , D) + I 3 (U n , D) 



where 



h{U n ,D) 



\U n 



xr ]n A(3{Qi) S(U n ) 



V \Qi\ 



h(U n ,D) = 



\Un\ 



jea, 



h(U n ,D) 

Here and below S(U) = S(U, D) for U C Z d and L> > 0. We have 



IQ^ES^Q^KviD), j e Z d , 
where w(-D) is the same as in (I16fl . By virtue of condition l|28[l it is clear that 

h{U n ,D) < \U n \~ 3 ES 2 (U n ) \Qi\ + W 2 t\s(Qj)s(u n )\ 

<i;(^){ixo(fo„)ii^r i +2|x (&n)r /2 ic/nr i/2 }^o as n ^^. 

For a fixed c > introduce the functions 

h\{x) — sign(ir) min{|x|, c}, h2(x) = x — hi(x), x G M. 
Given a nonempty finite set Q <Zl d let 

S(Q) = S(Q)/vlQ|. 

Note that 

2 



where 



h(U n ,D)< £ 4 P ' m) (t4,£) 



4^ m) (C/„,i?) = ^ A»,(3(Qj))M<S(Qj)) - Eh p {S{Q 3 ))h m {S{Q 3 )) 



For 6, n g N introduce the sets 



T^={seU n : inf ||s-i|| < 6} 

tedU„ 



(29) 



(30) 
(31) 



(32) 



where dll n = {j U n : 3q ^ U n such that \\j — q\\ = 1}. Put T n = Tn n \ n e N, where &„ meet 
condition (|28(l . Due to (|29|l one has 



4 1,2) ([/„,£>) + 4 24) (U n , D) + jf' 2) (C7 n , D) 
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<2\U n \- 1 E ( 2E\h 1 (S(Q j ))h a (S(Q j ))\ + EtZWQj)) 



:• 2( 2Ej/M.(,s-(/v I1 (^,))l^>(.V(/Cn(/v,)))i I Eh 2 2 (S(K (b n ))) + Sp^ig- 1 ^) 

1/2 



< 



4(v(D)E(S 2 (Ko(b„))l{\S(K (b n ))\ > c}) 
+ 2E(5 2 (^ (fon)l{|5(K (fc„))| > c}) 



where 1 is an indicator function. 
Condition (|29|l implies that 



+ 6|T„||C4r 1 t>(,D) 



iTnllt^l-^O as n^oo. 



(33) 



(34) 



It is easy to see that a family {S (Ko(b n ))}^—i is uniformly integrable. Consequently, taking into 
account llo-il) , for any e > we can find c = c(e) such that for all n large enough 



l£' 2) {U n ,D) + (U n ,D) + I {2 ' 2) (U n ,D) < e, 

furthermore, 

(I^ 1] {U ni D)f < \U n \- 2 E |cov^(5(Q,)),^(5(Q t )) 
In view of (|31(l we obtain the inequalities 



(35) 
(36) 



cov I h 2 (S(Q 3 )),h 2 (S(Q t )) 



\Un\- 2 E 

j,teU n ,\\j-t\\<4b n 

<2c 2 \U n \- 2 E e/i i(^(Qj)) 
j,tet7 n , 1||<46„ 

<2 M +V|ig- 1 |if (MK^)- 
Now condition l|l(Jfl with A > d entails the estimate 

\u n \- 2 E \<x"(%@(Qi)),ki@(Qt)))\ <c Q c 4 \u n \- 2 E \Q 3 \\Qt\\\j-t-2b ri 

j,t£U n ,\\j-t\\>4,b n j,t£U n ,\\j-t\\>4b n 



(37) 



I -A 



< c c A dA d \U n \- l \K {b n )\ 2 E ^ 1_A - 

r>26„ 

Taking into account (J3SJ) - (J2H1) , I® an d {El with A > 2d, we verify that 

h{Un, D) — > as n — > oo. 

Now observe that 

\Un\- 1 E iQj-r'E^w,) = ic/ n r l ic/„\T„iix (6„)r i E > s 2 (Xo(^)) 



(38) 



(39) 



+ic/„r i e iQ.r'E-s 2 ^). 

J6T„ 
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According to (|18fl the following relation is valid 

\K (b n )\^ES 2 (K (b n )) -» a 2 (D) as n -> oo. 
Due to and (J25J) we conclude that 

7 3 (C7 n , D) — * as n -> oo. (40) 
Relations l|3U)l. (JSUJ and JUJ yield lEU>- The proof of Theorem 2 is complete. 

5. Concluding Remarks 

Besides the concept of a functional reserve of an organ undergoing irradiation it seems desirable to 
use the models taking into account the geometrical configuration of survived FSUs (or/and cells). Here 
the concepts of random clusters appear naturally. In this regard we refer to a quite recent paper (cf. 
|17j). The stochastic models of disordered structures (involving point random fields) could be applied 
also to describe the damage volumes of irradiated organ. Note that it is possible to describe the relations 
between various FSUs (cells) in terms of random graphs and study the models where some vertices (or 
edges) are destroyed at random. An interesting problem is to find the optimal dose of irradiation taking 
into account not only the complication probabilities but the balance of conditions for irradiated organ 
(tissue) and its normal environment. Moreover, it is important to consider non uniform irradiation, 
another problem is to study a population of non- identical patients (see, e.g., [S3, El)- To conclude we 
mention a deep problem of constructing dynamical models describing the evolution of an irradiated organ 
in space and time. 
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